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SUMMARY 


This report deals with the problem of controlling the vibrations of large space 
structures by the use of distributed sensors and actuators. It specifically 
addresses the synthesis and implementation of Low -Authority Control (LAC) 
systems which arise when the above structural vibration controllers introduce 
only a moderate amount of damping (10 to 25 percent) in the structure. 

The LAC gain synthesis is achieved by an exact algebraic optimization of a 
weighted quadratic cost function based on the fundamental LAC root-shift 
prediction formula. The cost function is the sum of two disparate design 
objectives, making its optimization a special case of a pareto-optimal design [1], 
To broaden the applicability of the LAC synthesis procedure, the LAC theory 
has been generalized to arbitrary linear time-invariant systems which include 
general filter equations. The ensuing LAC synthesis procedure is then 
applicable to structures which already have some natural damping or some 
embedded attitude or structural control systems. In turn, this makes it possible 
to formulate and implement a sequential LAC synthesis procedure in which the 
perturbation method is applied sequentially in a manner analogous to the Newton- 
Raphson method. In this sequential procedure, the small amount of control 
synthesized at each step is embedded in the total dynamics , and the total 
system complex eigenvectors are recomputed for the next step. 

The robustness of LAC active damping systems is examined and analyzed for 
realistic (finite-bandwidth) actuation systems. The stability characteristics of 
idealized and increasingly more realistic "structure -with -actuators” models is 
analyzed, leading to the fundamental LAC Stability Theorem for already damped 
structures equipped with both (finite-bandwidth) active and parallel-mounted 
passive actuators. 
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The LAC synthesis procedure is applied to the design of LAC controllers for 
vibration control of simply- supported rectangular plates for which frequencies 
and mode shapes are given analytically. Two aspect ratios are considered: 

(1) a 3:1 ratio plate, for which it is shown that at least 10 percent damping 
can be achieved in the first five modes with only three colocated sensor /actuator 
pairs, and (2) a ^5/3:1 ratio plate, for which a similar objective is met with 
five sensor /actuator pairs. The special difficulty encountered for the second 
plate is the double-root frequencies of the fourth and fifth modes, leading to 
eigenvector indeterminacy. It is shown how the LAC sequential synthesis 
resolves that difficulty. 

A newly conceived linear actuation device, the pivoted proof-mass (PPM) damper, 
is described and analyzed. It consists of a flexure-pivoted mass (proof-mass) 
driven electrodynamically in an ac-coupled fashion, dissipating structural 
energy through electrical network heat losses. Detailed scaling laws are pre- 
sented for feasibility of extrapolation to large space structures. 
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INTRODUCTION AND OVERVIEW 


In an earlier companion study [2] , a general theory was developed for a class 
of control systems called Low-Authority Control (LAC) systems. When applied 
to structures, these are structural vibration control systems consisting of 
distributed sensors and actuators with limited damping authority, i.e., the 
control system is allowed to modify only moderately the natural modes and 
frequencies of the structure. This basic assumption, combined with Jacobi's 
root perturbation formula, leads to a fundamental LAC formula for predicting 
algebraically the root shifts produced by introducing an LAC control system. 
Specifically, for an undamped, open-loop structure, the predicted root shift 
dA. n is given by the quadratic form 

dA ^ - i ^ D 0 0 

n 2 ^ ar v an ^rn 

a,r 

where the coefficient matrix D ar is a matrix of (damping) gains, and 0 ^, 0 rn 

th 

denote, respectively, the values of the n mode shape at actuator station a 
sensor station r. 

Formula (i) may also be looked upon as (a set of) equations for the unknown 
gains D ar if the dAj^ are considered to be known (desired) root shifts or, 
equivalently, desired modal dampings. While an exact "inversion" of (i) does 
not generally exist, least-squares type solutions can be devised to determine 
the actuator control gains D &r necessary to produce the required modal damping 
ratios. This determination of the gains is the synthesis of LAC systems, and is 
the basic problem addressed in the present study. In order to make such an 
LAC synthesis as useful as possible for practical application to large space struc- 
tures, the LAC "scenario" has been generalized to arbitrary linear time-invariant 
systems of the form given in Eq. (1), p. 13, which, in this application , 
corresponds to already damped (controlled) structures equipped with 
distributed actuator /sensor systems, and for which the basic plant also 


(i) 


and 
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includes filter equations which may be used to model either sensor dynamics, 
actuator dynamics , state estimators , or any combination of them . In this most 
general case, the analog of formula (i) is the bilinear form 

_ 1 V .R 

dA.=-rr>D 0 0 

n 2 arn an rn 

a,r 

where the terms D = D (n) now depend on the mode number n and are 

arn ar ^ tJi 

generally complex, and where » 0^ n are, respectively, the a and r com- 
ponents of generalized actuator and sensor mode shapes 0^ and 0^. 

When the eigenvalues of the filter become large compared to those of the 

structure, it can be shown that the terms D become independent of n, and 

arn 

(ii) then reduces to 


d\ 


n 



a,r 



rn 


(iii) 


which is the most useful formula to use when sensor, actuator, or filter dynamics 
can be ignored or are already embedded in the system. Formula (iii) applies to 
already damped structures and is the basis for the general LAC synthesis pro- 
cedure (LACSYS) and its iterative version (sequential LACSYS) discussed in 
the first chapter of this study. 

The synthesis procedure is based on optimization of a weighted quadratic cost 
function related to (iii), and, in its implementation, is qualitatively similar to 
the classical optimal control gain synthesis procedures where a "dialogue" between 
the user and the free parameters of the system (in this case the weights in the 
cost function and the desired target modal damping ratios) is necessary to ensure 
that the control objectives are met. Even though synthesis of the control gains 
(for each set of specified weights and desired dampings) is achieved by exact 
linear algebraic processes, the resultant root shifts will only approximate the 
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desired ones because of the approximations inherent to the perturbation method 
employed. To improve the accuracy of the process, the perturbation method is 
applied sequentially in a manner analogous to the Newton-Raphson method. At 
each step, the small amount of damping due to the synthesized gains is "embedded" 
in the total dynamics, and the total system (complex) eigenvectors are recomputed. 
Because the synthesis procedure is applicable to already damped systems, it 
can indeed proceed in small steps such that the eigenvector shifts corresponding 
to the incremental damping added at each step remain small. This is particularly 
useful for the case of multiple roots for which the corresponding eigenvector 
indeterminacy produces large eigenvector shifts when damping is introduced. 

The second chapter deals systematically with the robustness of LAC active 
dampers . Robustness is that quality of a controller of remaining stable in 
the presence of parameter variations in the structure and/or control system 
parameters. For example, a stability condition independent of mode shapes and 
frequencies , such as the negativity of all root shifts in the fundamental LAC 
formula (i) for the case of colocation (a = r) , is an example of robustness 
when idealized (i.e., infinite bandwidth) actuators and sensors are used. 

The generalized LAC theory developed in the first chapter makes it possible 
then to examine various stability characteristics of idealized and increasingly 
more realistic (finite-bandwidth) systems, starting from ideal dampers in un- 
damped and damped structures , actual active or actual passive dampers in 
undamped and damped structures , to the final case of both actual active and 
passive dampers in damped structures. The different cases considered lead 
to the final Stability Theorem for LAC active dampers in the presence of both 
natural structural damping and parallel-mounted, device-implemented passive 
damping. 

The third chapter is devoted entirely to an application of the LAC theory in 
order to illustrate the LAC synthesis procedure. To this end, the structural 
model used is that of a simply-supported rectangular plate. For this model, 
the availability of analytic expressions for its frequencies and mode shapes 
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alleviates the need for finite-element models and their attendant numerical data- 
processing pitfalls. The 3:1 aspect ratio plate is treated first, and it is shown 
how three actuator /sensor pairs can easily produce (at least 10 percent damping 
in each of the first five modes. The sj5/3: 1 aspect ratio plate is treated next 
its characteristic feature being that its fourth and fifth modes correspond to a 
double-root. The sequential LAC synthesis procedure is applied again in this 
case and a complete tabulation is given for the desired dampings, modal weights, 
and cumulative LAC gains occurring at each step. For this plate, 5 actuator/ 
sensor pairs are used to produce (at least) 10 percent damping in each of the 
first five modes. 

The fourth and final chapter deals with the mechanization of LAC systems, and 
focuses in particular on the proof-mass actuator as the device closest to the 
linear analog of the gyrodamper studied previously in [1]. Proof-mass actuators 
operate by the principle of inertial reaction : they apply a force to a small mass 
(proof-mass) and the d'Alembert reaction force is transmitted back to the struc- 
ture. While this concept is quite straightforward, its practical implementation 
leads to numerous mechanical difficulties (e.g., friction, suspension stiffness, 
nonlinearities , etc . ) . 

Rather than attempting to model these various pathologies (which most often are 
not even reproducible in a systematic way), the approach was taken to place 
a control system around the actuator itself in order to force it to behave in the 
desired linear manner. Furthermore, a much cleaner mechanical design was 
obtained by approximating linear motions via a small-angle pivoting device, the 
pivoted proof-mass (PPM) actuator. An electrodynamically driven prototype 
was built and is described in this chapter, and detailed scaling laws are pre- 
sented for feasibility of extrapolation to large space structures. 


LIST OF SYMBOLS 


actuator index 
rectangular plate length 

filter dynamic matrix 
rectangular plate width 

PPM actuator composite center of mass coordinate 
unloaded PPM actuator center of mass coordinate 
filter input matrix 

active control gain of a generic damping actuator 
active control gain of damping actuator a 
active contribution to damping 
control gains matrix 

PPM actuator attachment point coordinate 

vector of selected damping gains for LAC synthesis 

i ** 1 component of d 

root shift of the n th root under LAC control 
damping gains matrix or scalar gain 
equivalent damping gain matrix with filtering 
element a,r of D 
element a,r of D 

n 

Young’s modulus 

actuator reaction force on structure 

maximum PPM actuator output force 

optimal PPM actuator output force 
PPM actuator electrodynamic force 
PPM actuator maximum electro dynamic force 



dynamic matrix 
input matrix 

rectangular plate thickness 
measurement matrix 
free index 

imaginary symbol y~l 

loaded PPM actuator centroidal inertia 
unloaded PPM actuator inertia w.r.t. pivot 
sectional inertia of uniform beam 


free 

index 


LAC 

cost function 


PPM 

center of mass 

scaling parameter 

PPM 

output force 

scaling parameter 

PPM 

inertia 

scaling parameter 

PPM 

mass 

scaling parameter 

PPM 

radius 

scaling parameter 


PPM actuator proof-mass coordinate 
PPM actuator characteristic length 

th 

left eigenvector corresponding to n eigenvalue 
uniform beam length 

modal index for rectangular plate 

mass of actuator proof-mass 
th 

i modal index for rectangular plate 

rectangular plate mass 
uniform beam mass 

passive control gain of generic damping actuator 

th 

passive control gain of a Ln damping actuator 
passive contribution to damping 





q 

% 


r 

r 

r 


R 

R 

max 


s 



ratio of roll-off time constants 

ratio of passive to active damper time constants 
sensor index 

radius of PPM actuator electro dynamic motor 
rectangular plate aspect ratio 

right eigenvector corresponding to n th eigenvalue 
ratio of active-to-passive damping 
maximum value of R ensuring robustness 

Laplace Transform variable 

matrix constructed with the mode shapes for LAC synthesis 
element of S 

PPM actuator torque output 


u 


input vector 


PPM actuator velocity 


W 

W 


R 


x 

x 



root shift weighting matrix for LAC synthesis 
modal weights 

abcissa of a point on the rectangular plate 

relative displacement of the actuator proof-mass 

normalized abcissa of a point on the rectangular plate 

t h 

normalized abcissa of a L actuator on the rectangular plate 
normalized abcissa of r sensor on the rectangular plate 
ordinate of a point on the rectangular plate 
measurement vector 

inertial displacement of a point of a structure 

normalized ordinate of a point on the rectangular plate 

th 

normalized ordinate of a LI1 actuator on the rectangular plate 

th 

normalized ordinate of r sensor on the rectangular plate 
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I 


z 

z 

z. 

1 

z 

mn 

Z 


( 5 .. 

On 


0an 

^rn 


filter state vector 

vertical displacement of a point of the rectangular plate 

th 

vertical displacement of the 1 rectangular plate mode 
vertical displacement of the (m,n) rectangular plate mode 

function of active and passive roll-off time constants. 

n tJl root of the beam frequency equation 

Kronecker symbol 

structure n 1 real mode shape vector 

th 

mode shape amplitude at a actuator 

th 

mode shape amplitude at r sensor 




complex ’’actuator mode shape’’ 


vector 


> R 

n 


*an 

R 

0rn 


n 


n 


© 


© 


m 


nth complex ’’sensor mode shape” vector 
component of 0^ at a** 1 actuator 

component of 0^ at r*^ sensor 

generic damping ratio 
active damping ratio 
damping ratio of n th mode 
passive damping ratio 

n** 1 eigenvalue, root 

proof-mass actuator angular travel 
maximum © 
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time-constant 


poles of active LAC damper 
pole of passive LAC damper 

surface mass-density of rectangular plate 
n** 1 frequency 




LOW-AUTHORITY CONTROL (LAC) FOR STRUCTURES 
Generalized Low Authority Control (LAC) Theory 


The LAC theory originally developed for uncontrolled and undamped structures 
Refs. [2] and [3]. may be generalized to any linear time-invariant system of the form 


X = FX + Gu ' 
y = HX 
z = Az + By 
u = Cz 


( 1 ) 


where X is the state vector [q q] 1 (q: modal amplitudes), u the control 
vector, y the measurement vector, z the filter state vector, and F, G, H, 

A , B , and C are constant coefficient matrices . In this formulation , the z 
equation (filter) may be used to model either sensor dynamics (in which case 
z is the actual sensor output) , actuator dynamics (in which case z is the 
actual force /torque output), any type of linear filtering including state estimators 
(in which case z is an estimate of the state X) , or any combination of the above. 

There are two reasons for extending the LAC theory to these more general 
systems: 


(1) The introduction of the general matrices F and G (representing 
the dynamics of the system) allows the handling of structures which 
already have some natural damping or some embedded attitude or 
structural control systems. As will be seen subsequently, such 
generalization of LAC will make it possible to develop a sequential 
algorithm for the control synthesis. 

(2) The introduction of a "filter” equation provides a convenient way to 
model sensor or actuator dynamics or controller roll-off, and has led 
to a very important result concerning the stability of real-life active 
LAC dampers. 
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It will be shown that when the general LAC formulas are applied to undamped, 
open-loop structures , they reduce to the simpler forms obtained directly in previous 
studies (Refs. [1] and [2]). 

In order to apply the perturbation theory, the matrix C (or B) in (1) is 
assumed to be ,T small ,T in some sense, so that the shift in the n 11 root of the open- 
loop system will be (shown to be) given by 


dA s 
n 


1 \ D 0 A 0 

2 ^ arn an y 


A. .R 
rn 


a,r 


th 

where D is the (a,r) element of the matrix 
arn 


D =C(A-AI) 1 B 
n n 


( 2 ) 


( 3 ) 


and 0^ n and 0^ n are respectively the a** 1 and r** 1 components of the generalized 
actuator and sensor mode-shape vectors given by 


0 


A 

n 


T 
G L 


n 


(4) 



HR n 


(5) 


where L and R are, respectively, the left and right eigenvectors of F 
^ ^ thi 

corresponding to the n root A n of F. 

Since eigenvectors are only defined to within a (complex) scalar multiplicative 
constant, defintions (4) and (5) depend upon the choice of separate normalizations 
for L and R . However, Eq. (2) is valid for any normalization of these vectors 
as long as, together, they satisfy the condition R n = 2. This particular choice 
is made so that Eq. (2) be formally compatible with Eq. (16) of Ref . [2] and 
Eq. (27) of Ref. [3]. 
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II 


th 

The indices a,r in Eq. (2) refer, respectively, to the a input u controlling 

th ° 

the system of actuators and the r measurement y^ of the system of sensors.* 

In simple cases, the indices a and r may be directly associated with (respectively) 

single- degree-of- freedom actuators and sensors. In more general cases, however, 

a single input u may drive several actuators, and/or a single measurement y 
a r 

may be a linear combination of several sensor outputs. 

The derivation of Eq. (2) proceeds from first rewriting Eq. (1) as the single 
matrix equation 


r * 


r 



X 


F i GC 


X 

z 


i 

BH. A 


z 


( 6 ) 


and considering the open- loop (unperturbed) matrix 


and the perturbation 




d <A 


F 

“1 

O 

BH 

• 

A 


O 

i GC 

O 

i 

i O 


(7) 


( 8 ) 


In order to apply Jacobi’s formula (Refs, 
eigenvectors and of <A are needed. 


II] or 12] , Eq. (5)) the right and left 
As can be easily verified. 


*The suffix "s" for sensor is not used in order to. avoid confusion with the common 
usage of s in Laplace transforms. 
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these are 


K 



n 



where 

R' = - (A - X I) -1 BHR 
n n n 


Jacobi's formula is then written as 



(9) 


( 10 ) 


( 11 ) 


( 12 ) 


or 


dA = - L T GC (A - A I) 1 BH R (13) 

n n n n 

which leads to Eq. (2) when the definitions (3), (4), and (5) are used. 

It is interesting now to apply Eq. (2) to an undamped, uncontrolled structure, 
for the case where 

(i) Sensors and actuators are physically colocated* 

(ii) Sensors and actuators are consistent, i.e., are of corresponding 
types: translation / force , rotation /torque 

(iii) Sensors measure rates 


*This does not necessarily require that the sensor-to-actuator feedback is also 
colocated. 
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In that case, we have 


H = [$ ; 0] 

G T = [<£ | 0] 


( 14 ) 


where $ is the matrix of (real) structural mode shapes 0 . In addition, if we 
normalize the eigenvectors L n and R r so that (compare to Eq. (XI) of 
Ref. [2] ) 


R = [0 0.. .0 1 0.. .0 ' 0 0. . .0. . . 1/ico 0...0] 

n » n J 

L = [0 0...0 1 0...0 j 0 0...0...SW 0...0] T 

n 1 n 


(15) 


(for which R n = 2, as required) , then Eqs. (4) and (5) reduce to 


0 A =0 
v an an 


0 R =0 
rn rn 


(16) 


so that Eq. (2) becomes the quadratic form 


d X 


n 



D 0 0 

arn an rn 


a,r 


(17) 


The above formula is very similar to Eq. (16) of Ref. 1 1] , but the filter 
[third equation in Eq. (1)] introduces two main differences: 

(1) The coefficients D arn = !* ar ( n ) now depend upon the mode number n 

(2) They are now complex quantities so that, in general, the dA f s are also 
complex, i.e. , first-order shift in frequencies may be obtained along 
with damping. 
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When the eigenvalues of the filter (i.e. , the eigenvalues of A) become large 
compared to those of the structure, Eq. (3) shows that D arn becomes independent 
of n and thus Eq. (17) converges to the old Eq. (16) of Ref. [1] and corresponds 
to the reduced equations 


Dynamics : X 
Sensors : y 

Controls : u 


= FX + Gu 
= HX 
= -Dy 


(18) 


The above equations, along with definitions (4) and (5) for generalized actuator 
and sensor modes, and the revised Eq. (17) for large filter eigenvalues (large 
filter bandwidth) , i.e. , 


dX = 

n 




ar 


a,r 


A R 
0 0 
an rn 


(19) 


are a convenient set to use when sensor, actuator, or filter dynamics can be 
ignored or are already embedded in the system. Indeed, in this latter case, 
using the definition (6) forc4, Eq. (1) can be rewritten as 


— — 


- ~ 


~ - 

X 

• • • 

= <A 

X 

+ 

G 

z 


z 


O _ 


u 


= [O : l] 


u = Cy 


X 


( 20 ) 


which has the same form as Eq. (18). Equations (18) and (19) are the basis for 
the general LAC synthesis (LACSYS) procedure and its iterative verson (sequential 
LAC SYS), which are discussed in the next section. 
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Low Auth ori ty Control Synthesis (LACSYS)* 

Equation (19) of the previous section is the basic LAC prediction formula for 
the root shifts (hence closed-loop damping) produced by sufficiently "small" con- 
trol (damping) gains D . To synthesize the gains D , let (dx ) denote 

ftr 8.1* ri p 

the predicted root shifts given by the above formula, i.e. , 


(dx ) = - \ S D 0 A 0 R 

n 2 Zj ar v an rn 

^ a,r 

= function (D ’s) 
ar 


(19a) 


and let (dX^).^ (given numbers) denote the desired root shifts imposed by the 
LAC controller design. Then the gains D a r are chosen so as to minimize the 
weighted quadratic cost function 


J(D) 


y W [ ( dX ) (dx ) l 2 + y D 2 

^ n I n p n D ar 

n L a,r 


( 20 ) 


in which the modal weights W n help specify pole locations, and the term 
improves robustness** of the controller, as will be shown later. Since the cost 
function J(D) is quadratic in (dA^)^ and hence, because of Eq. (19a), also 
quadratic in the elements of D, the gains can be obtained algebraically by 
solving the linear equations 


for (the elements of) D. 


" 9J (D) _ Q " 
8D 


( 21 ) 


*The material in this section is a continuation and elaboration of the short section 
on Damper Design and System Robustness appearing on p. 21 of Ref. [1] . 

**Robustness is that quality of a controller of remaining stable in the presence of 
parameter variations in the structure and/or control system parameters. For 
instance, a stability condition independent of mode shapes and frequencies, such 
as the negativity of all root shifts in the fundamental LAC formula, is an expression 
of robustness. 


19 


To carry out this procedure, the double-indexed gain matrix D is relabeled 

ar 

as a singly indexed gain vector d. and the multiply indexed (generalized) 

A R 1 

modal coefficients 0 0 are relabeled as a coefficient matrix S . , using 

an rn „ m & 

the same correspondence between the single index i and the pair of indices a,r, 
used above. The single-index relabeling scheme (a,r)— i , for actuator and 
sensor labels which excludes those pairs for which D ar is chosen to be zero*, 
produces a correspondence between a,r, and i such that: 


d. 

1 


= D 


a.r. 
1 1 


m 


- - R 

7T 0 0 

2 'a.n y r.n 


1 A 

3u: 

1 


( 22 ) 


Finally, let W denote a diagonal weighting matrix whose diagonal elements are W , 
and for simplicity, denote the desired root shifts (dA. by - ^ n w n » considered 
as components of a vector - (£o>) . 

With the above relabeling, the weighted quadratic cost function in Eq. (20) 
can be rewritten as 


J(d) 



or in matrix form 


J (d) = [Sd - (£w)] T W [Sd - (£w)] + d T d 


(23) 


(23a) 


From Eq. (23) we can calculate the partial derivatives 3j/3d. (where j is an index 
having the same range as i) as follows: 


1 9J_ 

2 3d. 

3 


■ I»n 


I 


S .d. - £ w 
mi n n 


n 


n J 


+ d. 
3 


( 24 ) 


*For instance , when colocated feedback is desired , 


D =0 for a ^ r . 
ar 
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or equivalently in matrix form (I: Identity matrix) 


= (S T WS + I) d - S T W (M (25) 

Thus, the condition 8J/3d = 0 leads to the LAC gain synthesis formula 

(26) 

In general, W is a diagonal matrix. If we assume further that the weights W n 
are all equal, W becomes a scalar matrix, say W = wl, and Eq. (26) can then be 
written as 

d = (s T S + | i) 1 S T (fo>) 

so that 

lim d = (S T S)~ 1 S T (fw) 

w— °« (27) 

= S* (&>) 

t T -1 T 

where S = (SS) S is the least squares pseudoinverse of S. This is the 

direct generalization of Eq. (18) of Ref. [2]. 

v 2 

The purpose and effect of the term 2 D in the cost function, Eq. (20) , 

d)i ui 

merits further clarification. Its genesis was the purely empirical observation that, 

for colocated feedback, when an exact inverse solution exists (i.e. , when the 

number of actuator /sensor pairs is equal to the number of modes to be controlled), 

an arbitrary specification of the desired damping ratios almost invariably leads 

to grains D . with large magnitude and mixed signs. In other words, the result - 
aa ■ 

ing control system removes energy from some parts of the structure while adding 
energy to other parts in order to achieve the desired (specified) modal 

21 



I 



damping distribution. This situation is, of course, very undesirable for achiev- 
ing robustness because all dampings gains must be positive to ensure uncondition- 
al stability. 

It was further observed that some overall damping is always obtained by 
setting all the gains to some positive value of much lower magnitude, even though 
the resulting modal damping distribution is then quite different from the desired 
one. It appeared, therefore, that adding a penalty cost for large gains to the 
cost function J would tend to prevent the occurrence of negative gains and thus 
ensure robustness. This tendency can be seen directly in Eq. (26) and its 
limiting form, Eq. (27), for large weights. While, for large weights, nothing can be 
readily said concerning the sign distribution of the gains (components of the 
vector d) , for small (positive) weights W such that W = wl we have 


(S WS + I) 


-1 


= I - 


T 

wS S 


1 2 
+ 2 w 


T 2 
(S X S) - 


so that Eq. (26) reduces, to first-order in w, to 


T 

d s wS ( ) 


Since, for all n, £ u >0,W >0 and the elements of (i.e. , 

= 1/2 0^ R j are also positive, it follows that all the gains (components of d) 
are also positive, and this is precisely the robustness condition which the second 
term in the cost function J in Eq. (20) or (23a) is intended to bring about in the 
LAC gain synthesis. 


Remarks: In the above synthesis , the feedback (damping) gains D were chosen 

ar 

to minimize the cost- function J in Eq. (20) which is the sum of two disparate design 
objectives : the magnitude of the errors between predicted and desired closed-loop 
damping, and the magnitude of the feedback gains. In this sense, an LAC system 
is a special case of a pareto-optimal design [1] in which multiple-objective cost 
functionals are optimized : the design is pareto-optimal if none of the multiple 
objectives can be further improved without degrading any other. While there exist 
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numerical techniques for multiple-objective optimization by linear programming 
methods using equality and inequality constraints. Ref. [4] , the exact (and relatively 
simple) algebraic solution given by Eq. (26) is, for the case of LAC systems, a far 
superior synthesis method than the more general approximation techniques. 

Sequential Low- Authority Control Synthesis Procedures • 

The low-authority controller design synthesis method discussed in the previous 
section is valid for a general class of linear systems, which includes, in particular, 
structures in which damping may be present. These systems are described in 
state-space form by the matrix Eqs. (18), which result in the closed-loop dynamics 

X = (F + GDH)X (28) 

For sufficiently ,T small" controls (or damping) D , written now as 6D , let 

dF = G(<5D)H (29) 

be considered as a perturbation of the dynamics matrix F consistent with the 
assumptions implicit in Jacobi's formula (ef. Eqs. (4), (5), taf Ref. [2], The closed- 
loop dynamics Eq. (28) can then be written 

X = (F + dF)X (30) 

The perturbation dF of the dynamics matrix introduces root shifts dA of the 

it. 

n n (complex) root of F. When the corresponding eigenvector shifts are 

small, LAC theory establishes linear relationships Eq. (19) between the dA T s and 

<5D. Conversely, synthesis of the control gains <5D corresponding to desired values 

of the dA’s (or tu/s) is also obtained by linear algebraic processes shown in 

Eq. (26) where the d's are obtained from a single-index relabeling of the doubly 

indexed quantities D (now designated as small quantities 6D ). However, 

&r &r 
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because of the approximations inherent to the perturbation method employed, the 
resultant root shifts will only approximate the desired ones. 

To improve the accuracy of this process, the perturbation method is applied 
sequentially in a manner analogous to the Newton-Raphson method. At each step i, 
the small amount of control due to 6D.. is "embedded” in the total dynamics, and 
the total system (complex) eigenvectors are recomputed. The corresponding 
control matrix gain increment <5D^ then reduces the difference between predicted 
and actual closed-loop poles. This is possible because: (1) the linearity-implied 

additivity of the process allows one to rewrite Eq. (30) as 

X = [F + dF 1 + dF„ +. . .+ dF ]X (31) 

12 n 


where 


dF. = G(SD.)H 

and where the corresponding final gain matrix is given by 

D = <5D 1 + <5D 2 +. . .+ 6D n (32) 

after an n-step iteration; and Eq. (2) , because of the applicability of the LAC 
procedure to general systems , in particular to already damped structures . A 
geometric interpretation of the sequential LACSYS procedure is shown in Fig. 1. 

The sequential procedure is particularly useful for systems in which multiple 
roots occur. For such systems, the introduction of small amounts of damping 
generally results in large eigenvector shifts for which the LAC theory is no longer 
applicable. However, once a small amount of active damping has been introduced, 
the eigenvector indeterminacy (resulting from multiple roots) is removed, and 
further iterations will then make it possible to accurately synthesize the control 
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CLOSED-LOOP SYSTEM: 

X 

= [F.dF 1 + dF,...-dF n ]x 

WHERE: 

dF 

= G <5D. H 
1 

GAIN MATRIX IS: 

D 

~ 5D, + <5D, + ... <5 D 
12 n 


PREDICTED 



Fig. 1 Sequential Low- Authority Control Synthesis (LACSYS) 

gains. This method will be applied in the sequel to illustrate an LAC design for a 
simply-supported rectangular plate of aspect ratio such that double-root modes 
occur early in the sequence of modes ordered by increasing vibration frequencies. 


25 











ROBUSTNESS* OF LAC ACTIVE DAMPERS 


The generalized LAC theory makes it possible to examine various stability 
characteristics of idealized and incresingly more realistic systems (i.e. , systems 
with finite bandwidth) . The different cases considered lead to the final Stability 
Theorem for LAC active dampers in the presence of both natural structural damping 
and implemented passive damping. 

Ideal Dampers in Undamped Structures 

This case is only mentioned for completeness , since it has been largely 
documented in the original LAC papers Refs. [2] and [31. Briefly summarized, 
for a colocated and consistent** set of actuators and sensors using rate feed- 
back, the root shifts 


d\ = — ^-£/D 0 2 (33) 

n nn 2 ^ a an 

a 

are always toward the left of the ico-axis for all the modes if all the gains D are 

a 

positive. This robustness result obviously assumed that both sensors and actuators 
have infinite bandwidth, and also that the structure was initially undamped. Several 
departures from this idealization occur in actual practical cases and will be examined 
next. 


*See footnote, p. 19 

**See statement (ii) , p . 16 
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Ideal Dampers in Damped Structures 


Since some small amount of natural damping is usually found in any structure, 
or could have been already introduced by passive means , it is of interest to determine 
whether the addition of active dampers can still meet robustness conditions. In the 
case of an already damped structure for which £ Qn denotes the existing small 
damping ratio of the n** 1 mode, the dynamics matrix F in Eqs. (18) now takes the 
form (see Eq. (10), Ref. [3]): 


f 'X 1 

\ 1 
N 

N l 

- 2 t v w 1 

on*p , 

V. 


N. 

V 

V 

1 

I , 

1 

0 

■ J 


with its eigenvalues given by: 


(34) 


X = - t co ± iaj ' v/l - 
n on n n ¥ on 


(35) 


and its right and left eigenvectors are given by ^let co * - 1 - £ 


IT 

on 


R = 
n 


L = 
n 


r 

(0...1...0 j 0... - (e on « n + *■£)/«*• ••»>]' 

[°. . .i + i t on w n/v • - ° ! °- • • • - 0 ] 


(36) 
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T 

as can be verified by inspection. Note that L n R n = 2; see discussion below 
Eqs. (4) and (5). Using now Eq. (14), the generalized mode shapes Eqs. 
(4) and (5) for consistent sensor /actuator systems are given by 


0 R = 0 
rn rn 


an 


1 + i .— on - 0 , 




on 


an 


Thus, the real part of the root shift is 


i&KdX n ) s 


£ co = 

n n 


~ V D 0 0 

2 / . ar an rn 


a,r 


(37) 


(38) 


which remains always negative for colocated systems. It is to be noted, however, 
that in this case, there is an additional first-order shift in frequency due to the 
presence of the £ Qn term. 


Actual Passive Dampers in Undamped Structures 


Because of the flexibility inherent in all materials, a passive damper must be 
modeled at least as an ideal damper in series with a spring, as shown in Fig. 2.* As 
a result , the force /displacement relationship , instead of being a pure rate-feedback , 



Fig. 2 Passive Damper Model 

^Indeed, for high frequencies, the "dashpot" part is basically locked, onlv the 
springiness of the device remains. 
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has a transfer function of the form 


f(s)/x(s) = - ps/(l + t q s) (39) 

where p is the damping gain, and t = p/k where k is the spring constant. 

This means that the actual device behaves like a pure spring rather than a damper 
at frequencies large compared to 1/2 irr . In terms of control theory, the dynamics 
of such devices can be represented by the first-order lag (or first-order "roll- 
off") 1/(1 + t s), placed at the output of the idealized device. Equivalently, 
it can be interpreted as a perfect rate-feedback system [the ps term in Eq. (39)] 
where the sensed rate has been "filtered" by the first-order lag filter 1/(1 + ^s). 
This second interpretation is convenient for applying the generalized LAC formulas 
to this example . Indeed , the total system equations in this case are 


X = FX + Gu 

y = HX 

1 ► 
z = - — (z - y) 

o 

u = -pz 

Comparisons of Eqs. (40) with Eqs. (1) yields the corresponding Eq. (3), 
written for p = p: 


D 

aan 


1 — ir w 
o n 

2 ,2 


+ T 


LO 


n 


(41) 


and thus, for a system of passive dampers (with colocated sensing) having the 
same time constant t , we obtain from Eq. (2) 


^e(dA n ) = 


£ to s 

pn n 


1 + T" U 
\ O 


I p 

n/ 


0 

a an 


(42) 


a 
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which again is always negative. However, in contradistinction to the ideal damper 
case, the root shift now depends on the modal frequency a? n , and the passive 
damping £ becomes very small as ru n becomes large. 

Actual Active Dampers in Undamped Structures 

This is the most complicated case, but also the most interesting in terms of 
practical implementation of LAC systems. There are quite a number of differences 
between actual active dampers and the passive damper model previously described. 
These stem from the usually complicated dynamics of the actuator itself (e.g. , 
nonlinearities). However, most of these unwanted characteristics can be sup- 
pressed by various methods (e.g., servo-loops around the actuator itself, 
electronic compensators, etc.), but one characteristic always remains, namely, the 
finiteness of the bandwidth. While it is reasonable to assume, for the purely 
passive damper, that the response f(s)/x(s) for all practical purposes is not 
zero at very high frequency, the same assumption cannot hold in the active case 
where it is necessary to introduce, at the very least, another "roll-off" filter 
in the model. 

The frequency response of such active systems will be described by the trans- 
fer function shown in Fig. 3: 

f (s)/x(s) = (1 - + T-g (1 + ^ (c > o) (43) 

where t < and where c is the active control gain. For frequencies smaller 
than 1/2 tt t * the system behaves like a damper, then more like a spring, and 
a final roll-off occurs after l/2 7r T 2 with the response going to zero with a -1 slope, 
as shown in Fig. 3. Unfortunately, this last roll-off may introduce enough phase 
shift to create instabilities if there are modes present in the roll-off region. 

This problem can be treated by introducing a filter, as in the previous case. 


31 



log f/x 



Fig. 3 Frequency Response of Active and Passive Dampers (Bode Plot) 


This time the filter and control equations are 



(44) 


and application of Eq. (3) yields a result analogous to Eq. (42) for a system of 
identical active dampers (c & = c) with colocated sensing 


<#e(dA ) =• 
n 


U) ~ - 


1 - T T U)‘ 

1 2 n 


cn n 


2 ( 1 + T l U n) ( 1+T 2“n)? 


C 0 
a an 


(45) 


This last expression for the active damping £ exhibits the frequency-dependent 
term (l - o^j which changes sign for > 1/ Thus, if there are any 

modal frequencies beyond 2 ir/\] t t , the corresponding modes will always be 
destabilized. 

This very disquieting result appears to put an end to any thought of achieving 
robustness or even stability in active damping of structures. But, however hopeless 
it may seem , it is not really serious because nature and engineering ingenuity come 
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to help in two ways. First, the above result is valid only for undamped structures, 
and will undoubtedly change if natural damping is taken into account. Second, as 
is graphically suggested by the curves of Fig. 3, a combined use of active and 
passive damping in the actuator itself could prevent this instability.* 

Actual Active and Actual Passive Dampers in Damped Structures 


This is the most general case, where a real structure, i.e. , one which has 
some small amount of natural damping, is controlled by a feedback loop whose 
transfer function is of the type (c > o, p > o): 


f (s) /x(s) 


cs ps 

(1 + T s) (1 + T 2 s) 1 + TS 


(46) 


where c and p are, respectively, the active (controlled) and passive gains of 
the damping actuator. For instance, c can be implemented via some electro- 
dynamic motor while the p may be realized by inherent back e.m.f. damping 
(produced by a low impedance power amplifier) or by an actual passive device 
mounted in parallel with the actuator. The time constants r and r are chosen 
to roll-off the active compensation well before additional poles due to actuator 
and amplifier dynamics are encountered. 

The total real part of the root shift ^e(dA^) = “ ^ n w n can now be obtained by 
summing up the contributions in the previous cases. For instance, assuming there 
are a number of similar actuators, with identical active and passive gains c = c 

d 

p =p , and with colocated sensing, we have:** 
a 


*This solution was originally proposed in 1978 by J. N. Aubrun and M. G. Lyons 
on the basis of root loci considerations. The ensuing stability theorem is due to 
Aubrun . 

**The case of nonsimilar actuators is conceptually straightforward, but computa- 
tionally quite involved since, for that case, (i = 0, 1, 2) must be replaced by 

T ia (actuator- dependent) and the corresponding terms can no longer be pulled 

out of the summation 2 in Eqs. (42) and (45). 

a 
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(47a) 


2 £ co 
n n 


2 £ co + V 
on n ^ 

a 


1 " T | C a 

i * n a 


j 1 +T o"n) + ( 1+T ^n)( 1+T 2“n). 


an 


where 


and 


2 £ co + 
on 


P T 


( X - T 1 T 2“n) 


n ( X ♦« ( 1+T K) ( 1+T 2“n) 


p f = Y p 0 2 = p Y 0^ = 2 £ co] 

Z-r a an ^ ^ 'an pn n | 

a a 

c’ = Y c 0 2 = c / 0 2 s 2 £ w 

/. a an Z-» an cn n 


(47b) 


(48) 


The question now arises: does there exist a relation between the parameters 
t q> Y t 2 » £ Qn > p > and c such that > when satisfied, 'the RHS of Eq. (47) is always 
positive for any co^? The answer is in the affirmative, and is given by an inequality 
which places an upper bound R max on the ratio R of the active /passive LAC 
damping ratios c T /[p T + 2 £ Qn w n 1 » and this in turn puts an upper bound on the 
achievable active damping ratio £ q . This inequality is characterized in the 
following theorem: 

LAC Stability Theorem (Aubrun): Unconditional stability of an LAC system is 

guaranteed if and only if, for each mode n, the active damping ratio £ cn i s less 
than a certain maximum £c n . This maximum active damping ratio in any of the 
modes within the bandwidth of an active LAC controller is proportional to the 
sum of the natural structural damping ^on anc * damping £p n introduced 

by a passive damper mounted in parallel with the actuator. Specifically, 


£ cn < £* = R a + £ ) 

cn max on pn 


where the value of the proportionality constant R is given by 

max 
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R maX s min (k + 2 -fk, r lV r o 2 ) 


where 


K = 


T l /T 2 + T 2 /T l + 2 


and where 1 /t 1/t^ are the poles of the active damper, and 1/r is the pole 

of the passive damper. 

Proof: Conditions are sought under which the RHS of Eq. (47b) is positive for 

all (robustness condition). Using the following definitions: 


<3 = V T 1 

% s T o /T l 

R -c«/[p' + 2^(1 + ,£«£)]► 


the robustness condition can be rewritten as 


(49) 


(1 + X 2 ) (1 + q 2 X 2 ) + R (1 - qX 2 ) (1 + q 2 X 2 ) > 0 
which simplifies to 

(q 2 - Rqq 2 ) X 4 - [r (q - q 2 )- q 2 - l] X 2 -fc R + 1 > 0 (50) 

2 

for all X. This is a quadratic (parabola) in the variable X , and a necessary 

4 

condition for positivity is the positivity of the coefficient of X 

r < q/q„ = T i V T o (51 > 
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We assume this condition can always be satisfied by the controller design. 
Inequality (50) will be satisfied if either of the following two conditions holds: 

2 

positivity of the coefficient of X , i.e. , 


R 


q + 1/q 

^ 9 

1 - % /( i 


or negativity of the discriminant of the quadratic, i.e., 


(52) 


a (q + 1/q) (1 - q /q)+ 2 

R 1 “ 2 2—5 — 

d+ %lq) 


R + 


(q + 1/q) ~ 4 

f f . 2 . . 2 

(i + q Q /q) 


< o 


(53) 


2 2 

To simplify the result, it will be assumed that q « q, i.e. r « t r In 

O O 1 z 

other words, the bandwidth of the passive damper is assumed to be much larger 
than that of the actuator. We define now 


Z = q + 1/q (54) 

2 

and since q Q /q <<: 1> conditions (52) and (53) become 

R < Z (52a) 

Z + 2 - 2 n/Z + 2 < R < Z + 2 + 2 \l r L + 2 (53a) 

The left inequality in Eq. (53a) may be ignored. Indeed, from Eq. (54) and 
Eq. (49) it follows that Z = (r^ +T 2 )/ T i T 2 ^ 2 and hence Z + 2- 2\/z + 2 <Z, 
so that if the left inequality in Eq. (55) were violated, it would simply be sub- 
sumed in condition (52a). Combining then inequalities (51) and (53a), the maxi- 
mum value of R is chosen to be 

R = min (z + 2 + 2 sJZ + 2, T.V 7 ' 2 ) (55) 

max \ ’ 1 2 o / 
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Now, within the controller bandwidth, we have w t. « 1, for i = 0,1,2. In 
this bandwidth, the robustness condition is thus achieved when 


c’ < R (p 1 + 2£ w ) 
max v on n 


(56) 


since R < R , and using Eqs. (48), we finally obtain 
max 


t < R a + £ ) 

en max on pn 


as was to be shown. 

Note. In the statement of the theorem, Z + 2 is denoted by K. The expression 

— 2 

K + 2 n/k" may also be written, after some manipulation, as (1 + 1/q) (1 +Jq) . 
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LAC ACTIVE DAMPING FOR SIMPLY-SUPPORTED RECTANGULAR PLATES 


Simply-Supported (SS) Rectangular Plates 

This chapter illustrates the synthesis of LAC active damping systems which 
was discussed at length in the last two sections of the first chapter. Because 
the strategy employed in any controller synthesis process is always dependent 
upon the characteristics of the structure to be controlled (i.e. , actuator /sensor 
selection depends on controllability /observability , which in turn depends on 
structural mode shapes) , it was suggested* that application of the LAC theory 
be made to simply- supported rectangular plates. The primary reason for this 
choice was that, for this class of structures and boundary conditions, there 
exist simple analytical expressions for both frequencies and mode shapes. This 
fact obviates the need for finite-element structural models and the ensuing 
numerical data processing which might obscure the controller synthesis process. 

A secondary reason for choosing simply- supported rectangular plates lies in the 
fact that, by proper choice of aspect ratio, double- roots can be made to occur 
early in the frequency-ordered sequence of modes. While there is, in fact, 
no unique eigenvector which corresponds to a double-root frequency but 
rather an eigenplane (invariant subspace of dimension 2) , plate theory provides 
in that case two linearly independent orthogonal eigenvectors which span the 
corresponding eigenplane. This makes it possible to examine the continuum of 
mode shapes corresponding to multiple roots and, if these occur early in the 
sequence of modes, to address the actuator /sensor location problem accordingly 
in a reasonably low-order model. 


Frequencies and Mode Shapes of SS Rectangular Plates 


Consider a simply- supported rectangular plate with the following parameters: 


Length: a (m) 

Width: b (m) 

Thickness: h (m) 

2 

Young’s Modulus: E (N/m ) 


Poisson Ratio: v 

o 

Surface Density: p(kg/m ) 

Flexural Rigidity: D (N-m) 

D = Eh 3 /[ 12 (1 - y 2 )] 


*By Dr. Larry D. Pinson, NASA Langley Research Center. 
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and let one of its corners coincide with the origin of a rectangular coordinate 
system (x,y,z) , as shown in Fig. 4, where a typical point of the deformed plate 

has coordinates x,y, and z (x,y). It 
is thus assumed that the deformation 



Fig. 4 Coordinates for Simply 
Supported Rectangular 
Plates 


field is one- dimensional , i.e. , all deforma- 
tions occur along the z-axis. 

It is well-known from plate vibration 
theory (e.g., Ref. [5]) that the vibration 
frequencies of such SS-SS-SS-SS rectangu- 
lar plates are given by the doubly-indexed 
family 



where r == a/b is defined as the aspect ratio of the plate. The plate deformations 
z mn which satisfy the boundary conditions and correspond to these frequencies 

are given by 

z mn ( x,y ) = a mn sin ( m7rx / a ) sin ( n 7iy/h) (58) 

where the amplitudes a mn are determined from the initial conditions of the problem. 

The two parameter families cu mn and z mn can always be linearly ordered by a 
single index i when the index pair (m,n) is ordered as (m. ,n.) (i = 1,2...) 
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according to some rule (e.g. , increasing frequencies). Using such a correspondence 
(m.,n^) — (i), we define 


CO. = GO I 

l m.n. I 
11 


(59) 


z i 


Z I 

m.n. 

l i J 


In order that the deformation functions (or surfaces) given by Eq. (58) be mode 
shapes, proper normalization is required. In particular, we must then have 
(i = 1,2...; j = 1,2...) 


If 

Plate 


z.(x,y) z.(x,y)p dx dy 


6 

i] 


(60) 


T 

which is equivalent to the standard normalization M4> = I" used m structures 

(see Ref. [6]) where $ is the matrix of eigenvectors (mode shapes), M is the 
generalized mass, and I is the identity matrix. Since the functions z^ are 
orthogonal, the normalization conditions of Eq. (60) reduce to 


b a 

P a^ n J J sin 2 (m-n sin 2 (nir g) dx dy = 1 (61) 

o o 

which leads to (for all indices m ,n) 


a 


ran 



2 

Vm 


(62) 
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where M is the mass of the plate. Using the transformation 


X = x/a 
Y = y/b 


the mode shapes for the plate can be written as: 


( 63 ) 


z (X,Y) = - sin m7rX sin n7rY 

mn Vm (64) 

x,y e [0,1] 

The affine transformation Eq. (63) maps the rectangular plate onto a square 
plate for any value of the aspect ratio. The surfaces z mn (X,Y), displayed in 
Fig. 5 for selected values of m , n thus represent (except for scale) the mode 
shapes of any rectangular plate . However , any specific ordering of the frequen- 
cies mn (and corresponding mode shapes z mn ) can only be obtained by specifying 
first an aspect ratio r in Eq. (57). 


LAC Active Damping For 3:1 Aspect Ratio Plate 


To illustrate LAC active damping for a plate for which (at least) 10 percent 
damping is required in the first five modes, with at most five colocated actuator/ 
sensor pairs, an aspect ratio of 3:1 is chosen. The physical parameters for this 
plate are as follows: 


a 

b 

h 

E 

P 

v 

D 

n/D/P 

M 


= 3 m 
= 1 m 

= 1 cm 

10 2 

= 6.9 x 10 N/m (aluminum) 

= 27.6 kg/m 2 for 1 cm thick 

= i/3 

= Eh 3 /[12 (1 - I, 2 )] = 6.46875 x 10 3 Nm 
= 15.3093 m 2 /s 

= pab = 82.8 kg 
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MODE 9 
120.24 Hz 



Fig. 5 (Concluded) 


and its first ten modal frequencies in Eq. (57) for 


2 


r 


a 2 /b 2 = 9 


are given by 


w. = (7r 2 /a 2 )\/D/p 9 n 2 ) 


(65) 


where the correspondence (m. , n.)— (i) is obtained as follows 
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Table 1. FREQUENCIES FOR 3:1 ASPECT-RATIO PLATE 


Mode No. i 

m. 

l 

n. 

l 


iO^TT (Hz) 

1 

1 

1 

10 

26.7 

2 

2 

1 

13 

34.7 

3 

3 

1 

18 

48.1 

4 

4 

1 

25 

66.8 

5 

_5_ 

_ 1 _ 

34 

90.8 

6 

1 

2 

37 

98.9 

7 

2 

2 

40 

106.9 

— 8 

3 

2 

45 

120,2 Double Root 

— 9 

6 

1 

45 

120.2 — J 

10 

_4_ 

_ 2 _ 

52 

138.9 


The choice of actuator /sensor locations is based on the locations of maximum 
modal deflection amplitudes of the plate , as shown by the modal controllability 
surfaces of Fig. 6. In J. N. Aubrun's Low Authority Control Theory, these arise 
from the fundamental root shift formula (e.g. , Eq. (16), Ref. [1]) for undamped 
structures . 


dA - 
n 


-iX D 


0 0 

ar an rn 


a,r 


( 66 ) 


where the coefficients are the damping (control) gains, with the indices a,r 

denoting, respectively, actuator and sensor locations on the plate so that 


0 = z (X , Y ) 

an n a a 


0 = z (X , Y ) 

r rn n r r 
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Fig. 6 (Concluded) 

for each vibration mode n. For a single colocated actuator /sensor (a = r) , the 

2 

r.h.s. of Eq. (66) becomes D (0 ) , and for a given root shift, the gain 

2 2 u 
D is minimum when (0 ) = [z (X , Y )] is maximum. Thus, when the 

aa all il 3 3 

actuator location (X , Y ) is varied continuously over the plate, the cor- 

2 a 2 a 

respondence (0 ) z (X,Y) generates a surface whose peaks represent 

ail il 

optimum actuator locations for minimum gain D for each mode n. These 

33 
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controllability surfaces are thus the "squares " of the mode shapes, and are 
shown in Fig. 6. Figure 7 shows five actuators placed very close to the plate’s 
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# ACTUATORS INITIALLY USED, THEN KEPT. 


Fig. 7 Actuator Locations on 3:1 Aspect-Ratio SS-Plate 

centerline Y = 1/2. The choice of these actuator locations is heuristic, and 

represents a compromise between exact placement on controllability peaks and 

available grid locations on the plate when a computationally tractable mesh size is 

used for the plate, e.g., 121 possible locations for actuator stations, as shown 

in Fig. 7. In general, one may tabulate the values of (0 ) versus a and n, 

an 

and examine the ensuing matrix for its largest entries. In the present case, 
examination of Fig. 8 reveals directly which actuators are most (or equally most) 
efficient for the first five plate modes. 

This matrix suggests that actuators No. 2 and No. 4 may be superfluous for 
damping the first five modes and this is indeed confirmed by the numerical 
synthesis of the LAC gains obtained by Eq. (26) of the LACSYS procedure. 
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ACTUATOR NO 


UJ 
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X 



3 





X 

4 

X 



X 


5 

X 


X 


D 


Fig. 8 Actuator Effectiveness Matrix for 3:1 Aspect-Ratio Plate 


Table 2 gives results of this synthesis for all five actuators. 


Table 2 LAC SYNTHESIS FOR FIVE ACTUATORS 


Mode 

OPEN-LOOP 

LACSYS 

LAC PREDICTED 
CLOSED-LOOP 

FULL EIGENANALYSIS 
CLOSED-LOOP 

No. 

Freq. (Hz) 

Weights 

Freq. (Hz) 

Damping (%) 

Freq. (Hz) 

Damping (%) 

1 

26.7 

7.50 

26.9 

11.0 

27.7 



2 

34.7 

0.75 

35.0 

11.5 

34.7 


CONTROL 

3 

48.1 

7.50 

48.4 

10.9 

48.7 

11.27 

►DESIGN 

4 

66.8 

0.90 

67.1 

10.0 

68.2 



_ _5_ 

90.8 

0.75 

91.3 

10.0 

88.9 

11.92, 


6 

98.9 

0 

98.9 

1.1 

98.8 

1.13' 


7 

106.9 

0 

106.9 

1.4 

106.6 

1.43 


8 

120.2 

0 

120.2 

1.4 

116.3 

5.52 

^SPILLOVER 


120.2 

0 

120.2 

1.4 

120.2 




138.9 

0 

138.9 

1.8 



137.7 
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The gains obtained with the above are (let D = D ) : 

aa a 


D 1 = 

1,830.0 

N-s/m ^ 

D 2 = 

1.9 

N-s/m 

°3 = 

544.0 

N-s/m 

a 

II 

86.1 

N-s/m 

D 5 = 

238.0 

N-s/m v 


( 67 ) 


so that S Id I = 2,700 kg/s. The gains D 2 and D 4 thus contribute relatively 
very little to the control effort , and removing actuators No. 2 and No. 4 yields 
the following results for the final closed-loop eigenanalysis : 

Table 3 CLOSED-LOOP FREQUENCIES AND DAMPING FOR 
ACTUATORS NO. 1, 3, AND 5 


Mode 

No. 

Closed- Loop 
Freq. (Hz) 

Damping 

(%) 

1 

27.49 

10.51 

2 

34.99 

10.58 

3 

48.58 

12.02 

4 

68.04 

10.03 

5 

89.07 

12.05 


with control gains 


a 

i-* 

n 

1,798.0 

N-s/m 

D 3 = 

486.0 

N-s/m 

°5 = 

350.0 

N-s/m 


( 68 ) 
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for which a measure of the control effort is 2 | D I = 2,635 kg/s. 

fl cL 

Remarks : (A) The desired modal dampings in (£co) of the LACSYS Eq. (26) 

used in the above synthesis procedure were 11 percent for the first five modes, 

and 0 for the remaining ones. The intent was to achieve at least 10 percent 

damping in the first five modes only, so no attempt was made to control the other 

ones. (This was already evident in the choice* of actuator /sensor locations which 

were such that modes Nos. 6, 7, 8, and 10 were essentially unobservable and 

uncontrollable, i.e. , 0 = 0 ~ 0 since the centerline Y = 1/2 is a nodal line 

— — — — - an rn 

for all these modes.) This restriction to the first five modes is quite arbitrary, 
and any number of modes could have been selected providing sufficiently many 
actuators are used. Indeed, it has been previously shown (p. 21, Ref. [1]) that 
in order to specify the damping ratio of N structural modes, N colocated 
sensor /actuator damping units are required. The fact that it is possible to use 
fewer than five actuators for the SS plate cannot be generalized to arbitrary 
structures; whether or not one may use fewer than units in the general case 
depends on the mode shape configurations. 

(B) Modes Nos. 8 and 9 correspond to a double root for which, strictly speaking, 
there is no (unique) associated eigenvector but rather an "eigenplane" in modal 
space. That is, any (pointwise) linear combination of the surfaces Zg (X,Y) and 
z g (X,Y) is again a mode shape surface corresponding to the double root “8 = "9 
at 120.24 Hz, and as an example, Figs. 9(a) and 9(b) illustrate the surfaces 
Zg + Zg and Zg - Zg, both of which are also vibration modes at that same fre- 
quency. The special difficulty with multiple root modes is that the presence of 
an actuator on one of their controllability surface peaks cannot guarantee any kind 
of damping performance because the vibration pattern simply shifts around so that 
the controllability peak reemerges elsewhere, away from the actuator, as can be 
seen in Figs. 9(a) and 9(b). From a mathematical point of view, the LAC per- 
turbation theory no longer applies because a small root shift can produce a very 

*The five actuators were placed slightly "below" the plate’s horizontal centerline, 
by a distance equal to 1/10 of the plate’s width, as shown in Fig. 7. The idea 
was to create a very small spillover into modes Nos. 6, 7,8, and 10 without any 
attempt to control them. In practice, actuator-mounting errors and resulting 
spillover will always occur. 
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large eigenvector shift. The method devised to deal with multiple root modes is 
a consequence of the sequential LACSYS procedure and is discussed next . 




Fig. 9 Example of Alternate Mode Shapes Corresponding 
to Double Roots 


LAC Active Damping For 5/3:1 Aspect Ratio Plate 


In the previous 3:1 aspect ratio plate, modes No. 8 and 9 were double- 
root modes. The study of such modes (if they were to be included without 
omitting previous ones) would then require a relatively high-order model, e.g. , 
of order 20 for the first 10 modes. In order to avoid higher order models, a 
new aspect ratio was sought such that double roots would occur last in a five- 
mode model. Inspection of Eq. (57) shows that this can be achieved for an 
aspect ratio r = \T 5/3 for which we obtain the frequencies shown in Table 4. 
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Table 4. FREQUENCIES FORn/573:1 ASPECT-RATIO PLATE 


Mode No. i 

m. 

l 

n. 

l 

/ 2 5 2 \ 

( m i + 3 n i ) 

w./27r (Hz) 

1 

1 

1 

2. 666. . . 

7.125 

2 

2 

1 

5. 666. . . 

15.141 

3 

1 

2 

7. 666. . . 

20.485 

— 4 

3 

1 

10.666... 

28.501 ^-1 n .. D , 
| Double Root 

— 5 

_2_ 

2 _ 

10. 666. . . 

28.501 — ' 

6 

3 

2 

15. 666. . . 

41.861 

7 

1 

3 

16.000. . . 

42.752 

8 

4 

1 

17.666... 

47,205 


where the plate parameters are the same as before, except that now b = 2.324 m. 
The first three mode shapes and the corresponding modal controllability surfaces 
are shown in Fig. 10. The fourth and fifth mode shapes, given by the surfaces 

z^(X,Y) = 2 /s/m) sin (3/rX) sin (7rY) 

z c (X , Y) = 2/JM) sin (2?rX) sin (2 ttY) / 

o 

correspond to the double root above at 28.501 Hz. These two surfaces are 
orthogonal eigenvectors corresponding to the double eigenvalue =w . i.e. , 


z 4 (X,Y) z 5 (X,Y)p dX-dY = 0 

Unit square 

and hence any (pointwise) linear combination of these surfaces corresponds again to 
an eigenvector at the same frequency. Thus, the surface z^(X,Y) defined by 

z 0 (X,Y) = cosO z 4 (X,Y) + sin0z 5 (X,Y) (70) 
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CORRESPONDING MODAL 
CONTROLLABILITY SURFACES 




CJI 

CO 


Fig. 10 First Three Mode Shapes and Modal Controllability Surfaces for J5/3:1 Aspect-Ratio SS-Plate 


is also a mode shape (with the proper normalization) at that frequency. As d 
varies, the one-parameter family of surfaces generated by Eq. (70) sweeps out 
the eigenplane corresponding to the double-root W 4 W 5* Figure 11 shows a 
sampling of the surfaces |z 0 (X,Y)j for d varying between 0 and 7r/2 in 15 deg 
increments , and these correspond thus to modal controllability surfaces of the 
continuous family Eq. (70). The peaks of these surfaces are shown in Fig. 12 
(for the first five modes), together with five actuator /sensor locations chosen 
for this aspect-ratio plate. (See Remarks below.) 

The final result of synthesizing a low-authority control damping system for 
the -J 5 / 3:1 aspect ratio plate is given in Table 5. 


Table 5. PERFORMANCE OF LAC SYSTEM FOR 
nT573:1 ASPECT-RATIO PLATE 


Mode 

No. 

Open-Loop 
Freq. (Hz) 

(FULL EIGEN ANALYSIS) 

Closed- Loop 
Freq. (Hz) 

Damping 

(%) 

1 

7.125 

8.391 

37.73 

2 

15.141 

15.184 

35.14 

3 

20.485 

19.275 

14.30 

-~ 4 

28.501 

26.800 

10.43 

— 5 

28.501 

28.352 

10.93 

6 

41.861 

41.478 

5.03 

7 

42.752 

43.197 

6.44 

8 

47.205 

45.363 

6.17 


The synthesis was carried out sequentially, and the values of the weights W n , 
desired damping ratios £ , as well as the LAC-predicted (P) and actual (A) 
damping ratios occurring at each iteration of the sequential process are displayed 
in Table 6. The objective, as before, was to obtain at least 10 percent damping 
in the first five modes in a "balanced” fashion, i.e. , without imposing excessively 
high damping in any particular mode. 
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Fig. 11 Continuous Family of Modal Controllability Surfaces for^fsTii:! Aspect-Ratio SS-Plate 
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14 27 *0 53 66 79 92 105 118 131 144 157 


O ACTUATOR (1, 5) 

X CPi: CONTROLLABILITY PEAK 

OF MODE #i (1, . 5) 


Fig. 12 Actuator Locations and Modal Controllability Peaks for \/5/3: 1 Aspect- 
Ratio SS-Plate 
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Table 6 SEQUENTIAL LOW-AUTHORITY CONTROL SYNTHESIS FOR 

ASPECT-RATIO SS-PLATE 


PERCENT CLOSED-LOOP DAMPING 

(P: PREDICTED BY LAC THEORY; A: ACTUAL OBTAINED BY FULL EICENANALYSIS) 


Mode 

No. 

Open- Loop 
Frequency 

Iteration 1 

Iteration 2 

Iteration 3 

Iteration 4 

Iteration 5 

Iteration 6 

Iteration 7 

Iteration 8 

Iteration 9 

Iteration 10 

Iteration 11 

Iteration 12 

























(Hz) 

P 

A 

P 

A 

F 

A 

P 

A 

P 

A 

P 

A 

P 

A 

P 

A 

P 

A 

P 

A 

P 

\ A 

P 

A 

1 

7.1 

6.1 

6.1 

20.4 

21. 1 

22.3 

- 

23.6 

23.5 

32.0 

- 

39.1 

39.8 

35.8 

34.5 

43.1 

43.7 

42.6 

42.4 

31.9 

- 

33.3 

33.7 

37.5 

37.7 

2 

15.1 

4.7 

4.8 

13.2 

13.9 

13.2 

- 

13.7 

13.0 

14.1 

- 

16.6 

14.8 

21.0 

25.2 

28.3 

27.8 

32.4 

32.8 

30.5 

- 

33.8 

33.8 

35.2 

35.1 

3 

20.5 

4.6 

4.6 

9.6 

9.2 

10.7 

- 

10.9 

12.6 

10.5 

- 

10.0 

11.8 

16.0 

16.7 

16.2 

16.3 

15.3 

15.3 

15.0 

- 

14.3 

14.2 

14.3 

14.3 

4 

28.5 

4.8 

4.8 

9.2 

10.8 

9.4 

- 

9.6 

12.7 

10.3 

- 

9.8 

5.2 

9.0 

10.8 

11.4 

9.4 

11.4 

11. C 

11.8 

- 

10.3 

10.4 

10.5 

10.4 

5 

28.5 

4.8 

1.3 

5.6 

4.1 

6.4 

- 

7.0 

4.4 

7.0 

- 

9.3 

14.4 

11.4 

8.4 

10.8 

12.9 

11.0 

10.7 

8.1 

- 

9.3 

9.7 

13.9 

10.9 


J 1 

5 


10 


10 


10 


23 


23 


- 


- 


- 


30 


30 


- 



£9 

5 


10 


10 


10 


11 


11 








_ 




_ 


Desired Dampings 
£ in (fu) Term 




10 


10 


10 


















of LACSYS m 
Formula (20) 1 ' 


5 





11 


11 


















f 4 

5 


10 


10 


10 


11 


11 


11 


12 


12 




- 


- 



f 5 

5 


10 


10 


10 


11 


11 


11 


11 


11 


- 


11 


11 



*1 

1 


1 


0.5 


0.5 


1 


1 


0 


0 


0 


1 


1 


0 


Modal Weights 

























W n in LACSYS 
Formula (26) 

W 2 

1 


1 


1 


1 


1 


1 


0 


0 


0 


0 


0 


0 


(Dimensionless) 

W 3 

1 


1 


1 


1 


I 


1 


0 


0 


0 


0 


0 


0 



W 4 

1 


1 


1 


1.2 


1 


1 


1 


1 


1 


0 


0 


0 



W 5 

1 


1 


1 


1.2 


1 


2 


1 


1 


1 


0 


1 


1 



25D 1 

487 


624 


1,090 


1,190 


766 


600 


934 


1,090 


1,180 


1,140 


989 


989 



L5D 2 

195 


995 


387 


317 


688 


996 


818 


977 


1,070 


997 


1,230 


1,309 


Cumulative 
LAC Gains (2) 

«°3 

281 


544 


1,040 


1,170 


988 


870 


1,830 


1,780 


1,880 


1,780 


1,850 


1,886 


(N-s/m) 

26D 4 

233 


609 


212 


106 


436 


714 


343 


452 


324 


202 


285 


348 



26D, 

5 

698 


592 


221 


288 


486 


630 


265 


492 


451 


243 


169 


246 

























F 

inal Gains 


(1; Desired dampings S n in (fw) of LACSYS formula ( 26) have no effect when corresponding weights W n are set to zero. 

Cl) The gain indices mutch the actuator names (numbers) shown in Fig. 12. In that figure, an actuator located at X,Y carries the grid label X + Y -1, e.g., Actuator #1 
is located at grid point No. 30. 


Remarks 


1) The sequential synthesis process was carried out for the first five modes only 
of an eight-mode model. As shown in Table 5 , the "spillover" to modes 6, 7, and 
8 produces about 5 to 6 percent damping in these modes. This spillover is a 
strong function of the relative location of actuators and controllability peaks for 
the "unmodeled "modes. 

2) Controllability peaks represent an optimum actuator location for given modes , 
but damping in those modes will also occur for nearby (nonoptimum) locations 
for somewhat higher actuator gains. For example, as can be seen in Fig. 12, 
actuator No. 2 is located at CP 4 (controllability peak of mode 4) and is as near 
to CP2 as the finite-mesh discretization of the plate will allow. (For this plate, 

a discretization of 12 x 12 = 144 subplates was used, with 13 x 13 = 169 grid points, 
including the simply- sup ported boundaries, available for actuator placement.) 
Placing another actuator at CP2 would not only be redundant for CP 4, but would 
provide excessive control on the plate’s horizontal centerline where also actuator 
No. 5 is located at CPI. 

3) For double-root modes, controllability peaks lose their significance because of 
the shifting modal pattern (Fig. 11) . The choice of actuator locations shown in 
Fig. 12 came as a result of some preliminary numerical experimentation based on 
placing each of the actuators at a controllability peak. For that choice, sequential 
synthesis invariably produced very high damping ratios in mode 1 and very low 
damping in mode 5, which has a nodal line (zero controllability) passing through 
the actuator located at CP 4. On the basis of these numerical observations, the 
actuators at CP2 and CPS were relocated on each side of CP5, as shown in Fig. 12. 
With this arrangement, actuators No. 1, 2, and 3 control the double-root mode 4-5, 
with actuator No. 2 also controlling mode 2 in a slightly suboptimal fashion. 

4) The sequential synthesis process shown in Table 6 reveals that mode 5 is still 
the most difficult to control. At iteration step No. 9, 10 percent damping (at least) 
is achieved in all five modes. The next three iterations serve only to distribute 
the damping more equally between these modes , and in particular to reduce the 
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42.4 percent damping 1 in mode 1. At this point, lack of a more specific control 
objective begins to matter because the synthesis could proceed in any of several 
ways. LAC systems are not particularly appropriate for specific (relatively high) 
modal damping requirements in selected modes , and other techniques (high- 
authority control, HAC) ought to be used. It is precisely in the context of 
such other controllers that LAC systems are useful, i.e. , they provide broadband 
damping (in many modes) to absorb the destabilizing spillover of the other mode- 
specific controllers. 

5) The sequential choices for the modal weights W n and desired dampings 
at each iteration step in Table 6 reveal a process which is characteristic of all 
optimal control synthesis processes, i.e. , the heuristic adjustment of weights 
(occurring in the cost function) vis-a-vis the control gains obtained. (For 
LACSYS, the desired dampings, in a sense, also play the role of weights.) 


59 




PROOF-MASS ACTUATORS 


Proof-Mass Actuator Concept 

The application of forces to a structure in space can be achieved by three 
principal means: 

1) Intrastructural actuators , pushing one element of the structure against 
another (e.g., member^dampers) 

2) Inertial reaction actuators, which create a force on the structure by 
reaction against a separate mass element which may either be discarded, 
as in the case of jets, or kept connected to the structure, as in the so- 
called "proof-mass” actuators 

3) Environmental actuators (e.g., magnetic, solar, etc.) obtained by inter- 
action of the actuator with the environment 

Although any of these methods could theoretically be used for the control of 
structures, practical considerations will help suggest the most likely candidate for 
each application. Type 3 may be discarded first because the force levels involved 
would usually be too low. Type 1 is very straightforward and is useful in many 
cases. However, it has two important drawbacks: first, because of its intra- 

structural nature, it tends to push energy around in the structure, and second, 
in the case of large structures, "member dampers" for instance will have to be of 
such size that their own flexible characteristics will come into play and add to the 
complexity of the problem. As for type 2, the jet solution may not be acceptable 
for various reasons : excitation of higher modes , difficulty of implementation , 
contamination of sensitive surfaces, fuel storage, etc. 

Thus, the best actuator is a momentum exchange device. In a previous study 
Ref. [1] , the use of control-moment gyros (CMGs) was discussed and analyzed. 

These devices are capable of imparting very large torques to a structure and thereby 
absorbing significant amounts of vibrational energy. Although there does not 
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exist an (exact) linear momentum counterpart of the angular momentum storage 
devices, the ,T proof-mass” actuator comes closest to the linear analog of a gyro. 

The important characteristic to be remembered is that the force f applied to 
the structure is the opposite of the force acting on the mass m of the actuator, 
and thus this mass will accelerate with the acceleration - f/m. Therefore, unless 
the force is reversed, the mass will continue to travel. This seems an obvious 
point, but it is an important distinction from the CMG case, where the gimbal 
angle stays constant when no torque is generated. This means that proof-mass 
actuators cannot be used for rigid- body control; their main purpose is to control 
vibrations, since in this case they will produce a force with no DC content. In 
fact, in practical implementations, these actuators will always be AC coupled to 
remove any possible bias , and also the mass will have to be physically restrained 
by a weak spring so it would not drift. 

This being established , proof-mass actuators can apply a very significant 
amount of vibrational forces to a structure, and conversely, absorb significant 
amounts of vibrational energy with the proper feedback loops. They have the 
great advantage of easy implementation, since they are self-contained and can be 
attached almost anywhere to a structure. Their only drawback is that they are not 
currently available in space technology. 

Active Versus Passive Damping 

The distinction between active and passive damping may seem very academic 
at first. Indeed, let y be the displacement of a point in the structure, then 
local damping may be introduced by applying a force 

f = - Dy 


Now , whether D is due to the viscosity of a dashpot liquid , eddy currents , or 
an actuator-electronic compensator combination, makes no difference to the result. 
For space applications, however, it is easier to use electromechanical actuators 
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than hydraulic or friction systems which are not very reliable. Also, values of 
D can be adjusted very easily in the electronics so that values may be achieved 
that are not obtainable with mechanical systems. Thus, except for these few 
points just mentioned, there does not seem to be a real issue. However, the 
situation becomes quite different when one considers how the force f is generated , 
and for proof-mass actuators it makes a difference. 

Consider the simple case depicted in Fig. 13. The actuator produces a force 
f on the structure and a reaction force -f on the mass m. The displacement of 
the structure at the actuator location is y, and the relative displacement of the 
proof-mass with respect to the structure is x. The dynamic equation of the proof- 
mass is then 

m (x + y) = -f . (71) 



Fig. 13 Proof-Mass Actuator Acting on Portion of Structure 

We will consider two cases: 

1) Passive Damping 

The force f is generated by viscosity, eddy currents, or other passive 
means, and may be written as 


f = Dk 


(72) 
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Using Eqs. (71) and (72), the transfer function may be obtained 


f/y 


- mDs^ 
ms + D 


( 73 ) 


Equation (73) does not exactly represent a pure rate feedback; only when s » D/m 
is it approximately one 


f/y = - Ds 

Nqw assume that damping needs to be introduced at some frequency <x)/2ir . 
Letting s = ico in Eq. (72) leads to 


f/y 


TA 2 

mD c o 
D + imco 


Thus, the imaginary part (i.e. , the one which will actually contribute to damping) 
is 


(f/y) 


damping 


r 2 n 2 i 

m D co 

-p. 2 , 2 2 

D + m co 


ico 


- D f (D) ico 


(74) 


The next question is: how much damping is achievable? Equation (74) shows 
that the function D’(D) vanishes both at 0 and infinity, and has a local maxi- 
mum at D = mco, where the equivalent damping gain is DJ^^ = 1/2 mco. Since 
the equivalent feedback law around the frequency co has the form 


f = - D f y 


(75) 


the damping introduced by the passive device is limited. 
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2) Active Damping 


Can this limitation be removed by active damping? If active damping meant 
merely to replace the mechanical gain D in Eq. (72) by an electronic one, the 
answer would obviously be no. However, active damping means also that sensing 
may be done differently. In this case, it is possible to sense y directly, rather 
than the natural "sensing" of x occurring in a dashpot. 

Thus, the control law, Eq. (72), can be replaced by 


f = - Dy ' (76) 

and in this case, there is, at least theoretically, no limit on what D could be, in 
opposition to the passive case in Eq. (75). 

Linear Proof-Mass Actuators 


A first attempt to build a proof-mass actuator was made using a linear electro- 
magnetic actuator (the Ling shaker model 102) which normally is used to shake 
structures from a ground base (see Table 7 for specs). 

Table 7. LING 102 CHARACTERISTICS 


Weight: 1.1 kg 

Max. Force Output: 7 N 

Max. Travel: + 3 mm 


This actuator consists mainly of a powerful, heavy, permanent magnet 
and a cylindrical coil free to move in the magnet gap and attached to an output 
rod, as shown in Fig. 14. A bellows-type suspension system maintains the coil 
centered, while allowing for small axial motions. Because the magnet and its 
casing are heavy, this assembly is used as the proof-mass, while the output rod 
is attached to the structure being controlled. 
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PERMANENT MAGNET 



Fig. 14 Simplified Cross Section of Ling Shaker Model 102 

There are two problems to solve before actually using this type of actuator: 

1) The bellows suspension is quite stiff and tends to stick at low force levels, 
so that the (force output) /(current input) transfer function is nonlinear, 
both in frequency and amplitude 

2) Despite this stiffness, it is yet incapable of holding the weight of the magnet 
in the gravity field, and is especially weak in off-axis directions 

In order to solve the second problem, a cradle had to be built to hold the output 
rod in a quasi-rectilinear path (Fig. 15). This is not really good mechanically because 
the coil does not exactly translate and thus may come into contact with the magnet for 
large motions. As for the suspension stiffness and stiction, which is aggravated by 
the cradle system because of misalignments and approximate linear motions, nothing 
much could be done mechanically, short of redesigning the whole actuator and 
centering system. 

However, one interesting solution was to wash out the actuator pathologies by 
introducing a feedback loop around the actuator. Ideally, the force output f should 
be a linear function of the input current. One possibility is then to measure this 
force (say by load cells) and feed it back to the actuator after comparison with the 
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STRUCTURE 


Fig. 15 Pivoted Cradle Suspension for ling Shaker 


input command. It turns out, however, that this method is not easy to implement 
practically. Another solution was thus developed which consists in controlling 
not the force, but the velocity of the output. 

Velocity-Controlled Actuators 

In order to eliminate undesirable dynamic characteristics of the actuator, the 
velocity of the output rod of the actuator is controlled by the feedback loop shown 
in Fig. 16. 


POWER 

AMPLIFIER 



VELOCITY 

SENSOR 


Fig. 16 Actuator Feedback Loop 
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Let ra be the mass of the proof mass , k the stiffness of the suspension 
system, and x the relative displacement of the output rod. The rod is attached 
to some structure where displacement is denoted by y. Thus, the equations of 
the system , in the s-plane are 


2 

ms (x + y) + kx = - f + G (v - v) 
v = sx 
y = z(s)f 


(77) 


where z(s) is a transfer function characteristic of the structure. If the structure 
is very rigid, z(s) is very small in the frequency range of interest (e.g. , the 
actuator mounted on a test bench). Solving Eq. (77) for v gives the transfer 
function 


v/v c 


k t ms 


1 + _____ ___ __ 

Gs (1 + ms z(s)) 


(78) 


Thus, for G large enough, the influence of m, k, and z are washed out and v 
basically follows the commanded velocity v . The bandwidth of such systems is 
basically given by 


~ G/m (79) 

and the minimum achievable velocity is given by 

v = f o /G (80) 

where f is the friction force in the actuator system. Thus linearity in both 
frequency response and amplitude is greatly improved by the velocity feedback. 
(Typical examples will be seen subsequently in. Fig. 20.) 
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In the case of the Ling shaker, it is possible to sense the velocity by 
measuring directly the back emf of the driving coil with an impedance bridge, 
as shown in Fig. 17. The advantage of this method is that it does not require 
additional hardware. However, the adjustment of the bridge may be difficult 
to achieve because of temperature drifts, especially when the power dissipation 
is large over a period of time. 



Fig. 17 Velocity- Con trolled (VC) Actuator Electronics for Ling 
Shaker Model 102 

Velooity-eontrolled (VC) actuators are quite different from force actuators 
and thus must be used differently. For LAC implementation, the usual control 
law in its simplest form is 


f = - Dy (81) 

If the structural displacement y is relatively small compared to the displacement 
of the proof-mass*, then 


*This is normally the case because the structure’s generalized mass is usually 
larger than the mass of the proof-mass. 
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f 


/V 


mx 


> 


(82) 


and thus Eq. (81) may also be written as 



(83) 


and integrating once leads to 


v = — y + constant (84) 

m J 

This shows that the velocity of the actuator must be proportional to the 
structural displacement . Therefore, when using a VC actuator, one should 
feedback position and not rate as was the case for force actuators. This is an 
unusual, but very important point to be noted. Another way to look at it is to 
consider the electronic compensator shown in Fig. 18. 

The commanded force f goes through a gain -1/m and is then integrated 
before going to the VC actuator. The velocity of the proof-mass is then 
v = - f c /ras and the force output is f = - msv = f . Thus the integrating 
circuit has transformed the VC actuator into a force actuator. In practical 
applications , however , it is important to remember that proof-mass actuators 
must be ac coupled, and thus the integration in Fig. 18 must be replaced by a 
compensation of the type s/(s + a)(s + b) , where a and b are small compared 
to the frequency of the lowest mode to be controlled. 

Pivoted Proof-Mass (PPM) Actuators 

Linear actuators are plagued by three main disadvantages : considerable 
weight, stiff guiding/suspension systems with inherent friction, and weakness 
of that system in transverse loadings. K. Silveira developed a new actuator 
concept at Lockheed in 1975 for fast mirror actuation to eliminate these pathologies. 
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STRUCTURE 

(1) DIRECT USE OF VC ACTUATOR 

(2) USE OF VC ACTUATOR AS A FORCE ACTUATOR 

Fig. 18 Electronic Compensation for Velocity-Controlled (VC) Actuators 

This concept was modified and adapted by J. Aubrun to the present pivoted 
proof-mass (PPM) actuator. In this actuator the linear motion of the proof-mass 
is approximated by a small circle of arc about a pivot point realized by a flexure. 
This type of flex-pivot has three advantages: 1) it is very strong in trans- 

verse loadings, 2) it has no stiction, and 3) it is mechanically extremely 
accurate. The actuation is obtained by a light electrodynamic motor. Figure 19a 
shows a picture of an actual PPM prototype and Fig. 19b a schematic of the 
device. Velocity feedback is used also with this type of actuator. Figure 20 
shows the open- and closed-loop transfer functions of the actuator. 
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a. Prototype of PPM Actuator 



b. PPM Actuator Schematic 
Fig. 19 Pivoted Proof-Mass (PPM) Actuator 
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Fig. 20 Open- and Closed-Loop Transfer Functions of a Velocity- Controlled 
PPM Actuator (Output Velocity Versus Commanded Velocity) 


The dynamics of the PPM actuator are not as straightforward as for the 
linear type. A dynamic model of the actuator is shown in Fig. 21. 



Fig. 21 PPM Actuator Dynamics Model 
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Because of the pivoting, both reaction force and torque are produced on 
the structure. Application of D’Alembert's principle shows that 


f^ = - mb ({r cos e + 0 sin 6 ) 

• 2 

f = mb ( 6 cos 6-6 sin 6) 

T = (I + mb 2 ) 6 - mbd (0 cos 0 - 0 2 sin 6) 


(85) 


Since the angular displacement 6 is usually limited to a few degrees, the 
above equations may be conveniently linearized. The longitudinal force f^ 
may be ignored so that the force and torque applied to the structure by the 
actuator are given by 


f s 

mb 6 

(86) 

T s 

(I + mb 2 - mbd) 6 

(87) 


where (see Fig. 21) : 


rf 

6 = §— ? ( 88 ) 

I + mb 

Equation (87) shows that it is possible to have no torque transmitted to 
the structure by choosing the attachment point at the distance 

d Q = (I + mb 2 ) /mb (89) 

Varying d above or below this value will change the sign of the output torque 
as well as its magnitude. The mass of the PPM actuator and its CM may be 
adjusted by changing the position and mass of a lead piece situated at the distance 
£ from the pivot point. 
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Let b Q , m Q , and I be, respectively, the position of the CM, the mass, and 
the inertia with aspect to the pivot of the ’’unloaded" actuator, and the 
mass of the lead piece with which it is loaded. Then 


m 


rri + 
o 1 


in b + m « £ 
o o 1 

m 


( 90 ) 


I + mb' 


I + m 1 £ 
o 1 


and using Eqs. (86) and (88) shows that 


f/f 


e 


( 1 + mb /m- £ \ 
o o 1 \ r 

1 + I o/ m l* 2 / * 


( 91 ) 


PPM Scaling Laws 


We must consider two important parameters in the design and use of PPM 
actuators besides the mass/inertia properties mentioned previously: 

(1) f = maximum electrodynamic force 

0IIl 

(2) 0 = maximum angular displacement 

Since the angular displacement 9 depends upon the frequency, optimal 
choices for the design parameters will depend upon the frequency range of 
application. From the two conditions 



em 


( 92 ) 


Ml 


£ 0 

m 


( 93 ) 
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and Eqs. (86), (90) and (91), two new conditions are found: 


f I 


1+ I o/ m l ^ 2 / 


7 f 
£ em 


( 94 ) 


| f | < (1 + m o b Q /m 1 g') m 1 ^w 2 6 


m 


(95) 


when w/ 2 7 r is the frequency of the output force. 

Design regions may be derived from these conditions. The quantities £!v 
and co may be chosen as main parameters and regions of possible values of f 
plotted as functions of them. Typical plots were obtained for a prototype PPM 
actuator whose characteristic parameters are shown in Table 8. 


Table 8. CHARACTERISTIC PARAMETERS OF 
PROTOTYPE PPM ACTUATOR* 


Mass 



Geometry 


Electrodynamic 

I = 1.5 l(f 5 

o 

. 2 
kgm 

b 

o 

= 0.016 

m 


m Q = 0.088 

kg 

r 

= 0.021 

m 

f = IN 

em 



e 

m 

= 0.067 

rad 



Two cases were studied, for two different design values of the proof- 
mass m^: 0.08 and 0.16 kg, respectively. Results for design #1 are shown 
in Figs. 22 and 24, and for design #2, in Figs. 23 and 25. Figures 22 and 23 
show the output force f as a function of f/r. Depending on the lever arm, 
more or less force can be obtained. The zero point (obtained for a negative 
value of £) corresponds to the composite center of mass being at the pivot 
point . 


^Prototype actuator (design #1) delivered to Dr. G. C. Horner, Langley 
Research Center. 
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Fig. 22 PPM Actuator Design #1 - Performance Region (f, £) 







The fact that the curves in Figs. 22 and 23 do not go to infinity for £ /r = 0 
but have a maximum is due to the inertia of the unloaded actuator which absorbs 
part of the energy going to the proof-mass m^. The straight lines represent 
the limits, at various frequencies, due to the angular travel limitation. It is 
thus seen that at a given frequency, there is an optimum value for £ /r, 
corresponding to the intersection of the curve and the straight line. If this 
value of £ / r is chosen for a given frequency , then the actuator will perform 
properly at all higher frequencies. For instance, if actuator (design) #1 has 
to control vibrations at or above 5 Hz, a value close to 3 is optimal for £/r. 

In this case, design #1 will lead to the values: 

£ = 0.062 m 

b = 0.038 m 

m = 0.168 kg 

I = 8.2 x 10~ 5 kg m 2 

The effect of increasing the mass m^ is to push up the angular travel 
limit, as can be seen by comparing Figs. 22 and 23, thus allowing the actuator 
to work at lower frequencies with basically the same force output. Obviously 
the increase in mass can be traded off against an increase of £1 r, but at 
the expense of the force level. For instance, to work efficiently at 2 Hz and 
above, actuator (design) #2 should have a value of £/ r of about 5.3. 

Another set of curves is shown in Figs. 24 and 25 where f is plotted 

2 

against the parameter (w/2?r) . Again the two limits (maximum electro- 

dynamic force, maximum angular travel) define a usable region which is, 
this time, a function of £!r. The frequency scale is also shown on the top 
of the plot. These plots show that above a certain frequency, the maximum 
force output becomes constant. The intersection of the two limits corresponds 
to the optimal frequency , for a given value of £ /r . 


I optimal 
for 5 Hz 


(96) 
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FREQUENCY (Hz) 



Fig. 25 PPM Actuator Design #2 - Performance Region (f , w) 




The optimal lever arm (<£/r) can be obtained from Eqs. (94) and (95), i.e., 



(97) 


This expression, which defines the abscissa of the intersection of the curve 
defined by Eq. (94) and the straight line defined by Eq. (95), is only valid if 
this intersection is on the right of the peak of the curve . Otherwise , the peak 
itself is the maximum value for f, in which case A /r is given by 


U/r) 


m 


m b 
o o 

m^ r 



m. I 

1 o 

2 V 2 
m b 

o o 




(98) 


Corresponding to Eqs. (97) and (98), the values of f are, respectively, 
given by 



[ mb + r m 
o o 


(^/r) ] co 2 0 

opt J m 


(99) 


f 

m 


1 1 
2 


f 

em 


( 100 ) 


As was done in a previous study (Ref. [l])for gyrodampers, scaling laws 
may be derived by considering geometric scaling. That is, calling L a 
characteristic length* of the system, the following scalings are assumed: 


*The parameter L is assumed to be dimensionless by letting L = 1 correspond 
to the prototype actuator (design #1) described in Table 8, p. 76. 


82 



II I 



X 


m 1 /rn 
1 o 


(101a) 


(101b) 


Thus Eqs. (97) and (98) become 


.- 1/2 *- 1/2 


K K„ (Leo) 
r f 


- K t 9 
I m 



(103) 


When Leo becomes larger and larger, (iJ/r)^^ ^ ecreases » while 
constant, and there is a critical value of Leo beyond which the optimal 
value for ( £! y ) is given by Eq. (103) instead of Eq. (102). Correspondingly, 
the maximum force output of the actuator will be given by Eq. (100) instead 
of Eq. (99). These last two equations may be rewritten, in terms of the 
parameters defined by Eqs. (101), as 
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(104) 


2 4 / K r K m K f 

co f . = Q K K, (Leo) + ( - r - 1 - K K 

opt m m b V q_ m 


m 


\ 1/2 
j(Lo)) J 


X 1/2 (La,) 3 


a, 2 f 


* K r K f K b 


-1 


m 


2 (l+X KjK^K 2 ) 1 2 -1 


(Leo) 


(105) 


Equations (102) to (105) constitute the basis for PPM actuator scaling and are 
graphed in Fig. 26 for the values of the scaling parameters in Eq. (101a) corre- 
sponding to the prototype PPM actuator with X = m^fm Q = 0.9 . These 
curves may be used, for design purposes, in the following way: Assuming that 

a value of f of the force output is required for some application at a fre- 

2 

quency equal to cu/2 7r , the value of the parameter is computed and 

2 

plotted on the right vertical axis. Using either the f + curve (if w f is 

“ 2 opi 

below the critical point) or the f curve (if w f is above it) , one obtains 
the corresponding value of L on the horizontal axis, and hence the value 
of the geometric scaling parameter L . Then, for that value of L , one 
determines the lever arm (f /r) , again using either the ( f / r ) Q p t or the 

(f/r) m curve (the latter being a straight horizontal line). An example of 
the procedure will be given in the next section. 


PPM-Damped Structure Scaling Laws 

The scaling of the actuator will ultimately depend upon the structure to 
be controlled. Scaling laws may be derived simply from the rate feedback Eq. 
(76) or (81) 


f = -Dy 

Thus, the force amplitude at frequency o; n is 
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(106) 


f 

n 


Dy co 
n 


where y n is the vibration amplitude. According 1 to Eq. (33), the associated 
damping produced by one colocated actuator /sensor pair is 


£ 

n 


1 T^S 2 

2 w n D ^n 


Thus the actuator force for achieving the required damping ratio £ must be 


f a 
n 


2 co 2 y 
n £ 

0 2 n 
n 


(106a) 


This relation, in conjunction with Eqs. (99), (100), (101), and (102) will 
allow the sizing of the required actuators. Although it is not practical or even 
possible to define scaling laws for a general structure /actuator system, some 
insight may be gained by considering a simple structure and observing the 
evolution of the different parameters as mass and length are changed. For that 
purpose, a free-free uniform beam damped about one axis by one PPM actuator 
is considered. The beam mode shape at the tip and its frequency are given, 
respectively , by 




> 


co =( 3 2 M 1/2 L 3/2 (El ) 1/2 
n n s s s 


(107) 


where L and M are, respectively, the length and mass of the structure, /3 
s s ] 

the root of the beam equation, and El the classical stiffness parameter. 
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Using Eq. (106) leads to a condition on the actuator force output of 
the form 


^ max 
2 n 


K y n £ n) 


(El ) 
s 


-3 


(108) 


The 100 m graphite- epoxy beam described in Ref. [1] p. 46 will be used here as 

an example of application of these scaling laws. Table 9 gives characteristics of 

this beam from which one obtains the quantities L , M , and El . Table 10 

s s s 

gives the value of . The following problem is addressed: 

Assuming that 10 percent damping is desired in the first mode of the beam, and 
that the vibration amplitude of the tip is initially equal to 1 cm , design a PPM 
actuator from the results obtained with the prototype. 


Table 9 GRAPHITE EPOXY BEAM PARAMETERS 


Length 

Outside radius 
Wall Thickness 
Young's Modulus 
Density 
Mass 

Sectional Inertia 


L = 100 m 
s 

r = 10. 55 cm 
e = 2.275 mm 
E = 3.45 10 11 N/m 2 
P = 1607 kg/m 3 

M = 239.7 kg 

5 -6 4 

I = 8.125 10 m 
s 


Table 10 FREE-FREE BEAM EQUATION ROOTS 


n 


1 

4.7124 

2 

7.8540 

3 

10.9956 

4 

14. 1372 

5 

17.2788 
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First, by applying Eq. (108), it is found that a force f ~ 0.7 N is 
required. Since the frequency of the first mode is 0.405 Hz, 

= 4.53 

4 

This value is far below the critical point (~4.2 x 10 ) , shown in Fig. 26. Thus 
the set of curves shown in Fig. 27 which correspond to that lower range will be 
used.* These curves, (A/r) ^ and f Q j., have been parameterized with X = m^/m Q 

varying from 0.125 to 4. The value of 0.9 corresponds to the prototype PPM 
actuator. Using this value of X first, the value of Leo is found to be approx- 
imately 7.3, corresponding to L~2.9. The (j?/r) curve shows an optimal lever 
ratio of 13. If instead, one chooses X = 4, , smaller values for L and £fr 
are obtained. These two designs are compared in Table 11, obtained from the 
graphically determined values of L and if r and from Eqs. (101). 

Table 11 PPM ACTUATOR PARAMETERS FOR 100 m Gr-Ep BEAM 


Actuator 

Parameters 

X = 0.9 

X = 4 

L 

2.9 

2.3 

£ (r 

13.0 

7.5 

r 

0.061 m 

0. 048 m 

£ 

0.792 m 

0. 362 m 

m 

o 

2.146 kg 

1.071 kg 

m i 

1.931 kg 

4.283 kg 

m 

4.077 kg 

5. 354 kg 

f 

em 

8.41 N 

5.29 N 


♦Figure 27, for x = 0.9, represents the same curves (&/r) t 
and f as appear in Fig. 26, but drawn in a neighborhood 
of the^rigin for vastly reduced values of a/f and Leo, and 
correspondingly increased values of (£/r). In addition, these 
curves in Fig. 27 are now part of one-parameter families 
obtained by varying the parameter %. 
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Fig. 26 PPM Actuator Design Nomograph: Force and Lever-Arm Versus Geometric Scaling Parameters 








0° Fig. 27 PPM Actuator Design Nomograph for the 100-m Gr-Ep Beam (x = n^/ri^ as a Parameter) 



The performance of the system for higher modes can be obtained easily by 
Eq. (106). If it can be assumed that the modal amplitudes vary with 1/w n 
(this corresponds to the modal content of a static deflection), then Eq. (106) 
shows that the maximum achievable value for £ n is a constant. 

As can be seen, the second design, with X = 4, is a bit heavier but much 
smaller in size and requires a smaller electrodynamic force. It can be shown 
that for this particular design, and at this low end of the range, the 
minimum weight design is obtained for X = 1/2 , in which case the total mass 
of the actuator is 3.977 kg. However, with m , the corresponding 

linear dimensions become prohibitive. 
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CONCLUSION AND RECOMMENDATION 


The analytical framework has been established for the design and gain synthesis 
of Low -Authority Control systems for large space structure vibration suppression 
and management. Design implementation of LAC systems have been illustrated 
for rectangular plates equipped with infinite-bandwidth sen sing /actuation systems, 
and robustness criteria established for realistic (finite-bandwidth) actuators. 

A newly conceived electrodynamically driven actuator concept is given, together 
with an actuator control system (actuator velocity control loop) designed to 
minimize spurious mechanical effects which arise with any physical actuator 
containing some degree of structural flexibility. 

i 

The inclusion of actuator and/or sensor dynamics in these types of study can 
only be meaningful in the context of an experimental program in which laboratory 
brassboards , equipped with sensors and actuators , are used to demonstrate the 
mechanical implementation and performance of such LAC vibration control systems. 
It is recommended, therefore, that any future work in this area be strongly 
coupled to a laboratory experimental program in which the practical implementation 
issues are addressed directly. By doing so, it is hoped that a new impetus will 
be provided in the aerospace industry for the conception, design, and manu- 
facture of badly needed actuator devices for large space structure vibration 
control. 
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